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Abstract 

The linear stability analysis of an optically thin plasma where a 
general reaction proceeds, including chemical relaxation time effects, 
is carried out . A fifth order dispersion equation (instead of the fourth 
order one resulting when such effects are neglected) is obtained. The 
new mode with the corresponding instability criterion as well as the 
modifications of the previous four modes and the corresponding insta- 
bility criteria, are analyzed. Generally, a further stabilizing effect on 
the unstable modes and an increasing of the damping of stable modes 
appear because of the second viscosity generated by the chemical re- 
action. The results are applied to: (1) a collisionally ionized pure 
hydrogen plasma heated at a constant rate per unit mass and cooled 
by free-free transitions, ionization, and e — H collisional excitations; 
(2) a diffused gas with metallicity Z, photoionized and heated by a 
radiation field, and cooled by excitation of hydrogen and heavy metal 
lines. 



1 Introduction 

It has long been recognized that the knowledge of the thermochemical equi- 
librium and stability of plasmas is a crucial point for understanding the 
origin of astrophysical structures at very different scales ( Field 1965 ). In 
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fact, the onset of thermal ( or more generally thermochemical ) instabilities 
can directly induce the formation of inhomogeneities or can set in physical 
conditions appropriated for developing other kind of instabilities, say for in- 
stance, the gravitational one. As a direct agent in forming condensations it 
has been considered in studies of the solar atmosphere (Priest 1987; Reale, 
Serio & Peres 1994 ), the interstellar medium (Field Goldsmith, Habing & 
Field 1969; Hunter 1970; McKee & Ostriker 1977; Flannery and Press 1979; 
Cowie, McKee & Ostriker 1981;), gas at high latitude and in outer disks of 
galaxies ( Verschuur & Magnani 1994, Ferrara & Field 1994 ), Lya clouds 
(Sargent et al. 1980; Black 1981; Bond , Szalay & Silk 1988; Baron et al. 
1989; Murakami & Ikeuchi 1990; Pettini et al. 1990; Duncan, Vishniac & 
Ostriker 1991); quasars gas ( Krolik, McKee& Tarter 1981; Mathews 1986; 
Mathews and Doane 1990; Goncalves, Jatenco-Pereira & Opher 1993 ), and 
as an indirect promoter of the formation of stars ( Glassgold and Lager 1976; 
Oppenheimer 1977; Sabano and Kannari 1978; Ibanez 1981 ), globular clus- 
ters ( Fall and Rees 1985; Murray & Lin 1990a,b ) and galaxies ( Sunyaev 
and Zel'dovich 1972; Gurevich & Chernin 1975; Zel'dovich & Novikov 1983 

)• 

On the other hand, the generalization of Field's (1965 ) basic study on 
thermal instabihty for including effects of chemical reactions has been car- 
ried out by Goldsmith (1970), Defouw (1970), Yoneyama (1973), Flannery 
and Press 1979; Ibanez & Parravano (1983); Ibanez & Mendoza (1990), and 
Corbelli and Ferrara (1995 ). However, in these works it has not been taken 
into account the fact that when the temperature, density and the other ther- 
modynamical quantities change, the position of the chemical equilibrium also 
varies, i.e. in a reacting gas the fluctuations occur on states out of thermo- 
dynamic equilibrium because, generally, the chemical relaxation time is not 
short enough to follow the change of pressure and density. So, if the restoring 
of the chemical equilibrium after an initial compression, for instance, occurs 
relatively slowly, it can not follow the compression. But the irreversible ten- 
dency of the concentrations to reach the equilibrium values corresponding to 
the new value of pressure and density produces an increase of entropy with 
the corresponding energy dissipation. ( Landau and Lifshitz 1987). The 
present work is aimed at considering such effects in a general reacting gas. 
Additionally, two particular applications will be carried out: (1) to the pure 
hydrogen plasma model worked out by Ibanez & Parravano (1983), and (2) to 
the photoionized plasma with mctallicity Z studied by Corbelli and Ferrara 
1995 , hereinafter references IP and CF, respectively. 
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2 Basic Equations 



For a fluid where a chemical reaction of the form ^iCi = is proceeding, 
the equations of gas dynamics can be written in the form 



AR 



dT p dp 



| + V.(pv) = 0, 



dt 



^ + BRTf + T, i) - -V-(/^VT) = 
dt dt p 



p{p,T,0 = R 



(1) 

(2) 
(3) 
(4) 
(5) 



where ^ is the convective derivative and the remaining symbols are defined 
in Table 1. 

For reacting fluids is useful to rewrite equation (3) in the form 

d^ <l>{p,T,0 



dt 



where 



and 



r{p,T,0 



t{p,t,o ' 



)p,i 



<P{P, T, ^) = —Qj^ 



(6) 



(7) 



(8) 



see for instance Landau and Lifshitz (1987), Zel'dovich and Raizer (1966), 
Vicenti and Kruger (1975). When the reacting fluid is chemically stable > 
0, r becomes the chemical relaxation time, otherwise, r is the e — folding 
time for chemical instability. 

The chemical equilibrium is deflned by the condition 



0(p,T,O = O, 



(9) 
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therefore, the chemical equihbrium value of ^ ( denoted by becomes a 
function of the density and the temperature which generally are functions of 
time, i.e. = ^p{t),T{t)]. 

On the other hand, and as it is well known, in an hypothetical infinite 
medium, equations (l)-(5) admit a steady state (uniform ) solution corre- 
sponding to a complete thermodynamical equilibrium, i.e. v = 0, and p and 
T satisfying the equation 

£(p,T,O = 0. (10) 

Equation (10) defines the thermal equilibrium. 

The stability of the above equilibria is usually analyzed ( in the linear 
approximation ) looking for solutions of equations (l)-(5) in the form of 
plane waves, superimposed to the equilibrium values i.e. 

Sij = iJo + i^'it) exp(ik ■ r) , (11) 

where iIjq is the ip value at equilibrium and ip'{t) is the respective time de- 
pendence of the disturbance on the variable ip. Substituting equation (11) 
in equations (l)-(5) and neglecting nonlinear terms in ip' one obtains the set 
of linearized equations 

^ + ipok-v' = 0, (12) 
Po-^+ikp' = 0, (13) 
r-^ + e - C*pP' + C*tT' , (14) 

P' = PpP + PtT' + Pi:^' . (16) 

Hereinafter, the subindex q denoting equilibrium values will be dropped. 

The linearized equation (14) takes into account the fact that when p and 
T change, the position of the equilibrium also varies ( Landau and Lifshitz 
1987; Zel'dovich and Raizer 1966 ). For a detailed analysis of this basic 
aspect of the problem, see Vincenti and Kruger (1975). 

As usually, one may proceed further looking for solutions of the dynam- 
ical variables v'(t), p'{t) and T'{t) ~ exp(A/'t). However, from equation (14) 
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follows that a solution for ^' ~ exp(A/'t) exists if and only if the amplitude 
of the disturbances ^i, piand Tiare related by the equation 



where pi and Ti are the amplitude of the density and temperature fluctuation, 
respectively. Therefore, equations (12)-(16) are reduced to two algebraic 
equations for pi and Ti, i. e. 



[{Afr + l){Af^+Ppe)+C.pP^e)]pi + [{Afr+l)pTk^ + UP^e]Ti = , (18) 

[{J\fr+ l){-^J\f+Cp) H*p{BRTJ\f+C^)]pi + [{Nt+ 1) {ARU+Ct+ —) + . 
p^ p 

CMBRTAf + C^)\Ti = 0, (19) 
which have non trivial solution provided that the compatibility condition 

aoA^ + aiA/^ + a2A^ + a3iV^ + a4^+a5 = , (20) 

holds, where 

ao^r^ , ai = t(2 + -jT^*t) + r^ch , 

02 = (1 + T^*^^) + '^^(2^1 + ^3) + r'^c^k'^ , 

03 = c{ki + ks) + + T^dki - k2)k'' , (21) 

7 7 

04 = ^k^ + T^[(2 + ^.pppiu)k, - (2 + ^.TTpiu)k2 + (1 - ^)k3]e , 
7 7 ?*T 

«5 = ^[(1 + e*pP)"z^)A;i - (1 + UTHk2 + (1 - |^)A;3]A;2 . 
7 ?*r 

In previous works ( Goldsmith 1970; Defouw 1970; Yoneyama 1970; Flan- 
nery and Press 1979, IP and CF ) where the change in the chemical equilib- 
rium position was not taken into account, the characteristic equation found 
was a fourth order , instead of a fifth order polynomial, equation (20). Note 
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that the way of considering the disturbance in the chemical parameter ^ 
in eqs.(14) and (17) is different from that of the above papers. In fact, in 
those papers the hypothetical initial "equilibrium state" is considered frozen, 
as that in absence of relaxing processes and the chemical fluctuation is as- 
sumed to follow that of pressure,density and velocity. However, this is only 
an acceptable first approximation. In deed, once the initial equilibrium is 
disturbed, an irreversible process set in and the chemical evolution tends to 
a new equilibrium corresponding to the new values of pressure and density. 
Therefore, it is not matter how tiny this effect is, it changes qualitatively the 
problem. The effect of finite value of the relaxation time is equivalent to think 
on a second viscosity (Landau & Lifshitz 1987) and this is the last physical 
reason by which fluctuations with large wave numbers can be damped in a 
fluid where no other irreversible process, except the chemical reaction, goes 
into play. This aspect of the problem will be again considered later on. 

3 Instability Criteria 
3.1 Asymptotic Cases 

Equation (17) gives the correct amplitude for the chemical parameter in 
both asymptotic limits: no reacting gases, i.e. when r oo, —>■ 0, and 
in equilibrium flow for which r — > and — > C*pPi + C*tTi, i.e. when the 
chemical reaction is so fast that it instantaneously adjusts to the chemical 
equilibrium values and therefore, the changes of density and temperature 
occur in chemical equilibrium. Note that in the equilibrium flow limit, = 0, 
T — but {d^/dt)eq. 7^ 0, contrary to the simple equilibrium ( as in a closed 
system with fixed conditions) for which = and {d^/dt)eq. — 0, case 
considered in the mentioned references at the end of the previous Section. In 
practice however, and from the physical point of view, the above asymptotic 
limits correspond , respectively, to the case when the chemical time scale is 
much longer, ( \Nt |>> 1) and much more shorter ( \Nt |<< 1 ) , than the 
characteristic time of the density and temperature fiuctuation. 

On the other hand, as it can be readily verified, when r oo, the 
dispersion relation (20) reduces to the corresponding Field (1965) equation, 
and when r — > 0, equation (20) reduces to 

a2Af^ + asAf^ + aj^f + 05 = , (22) 
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where 02 



02 = 1 + ^T^^T ; 03 = c{ki + ka) ; 04 = ^-^/c^ , (23) 
A 7 

and 05 stands as in equation (21). 

By applying the Hurwitz criteria to the equation (22) one obtains that 
an instabihty sets in when any of the following relations is fulfilled, 

>CT + e*TA + — <0, (24) 

T T T nk'^ 

(i+r/xi/e*T)[(A+sre*T)>Cp+— (i-Sp/xe*p)>Cr+(Ae*p+— e*T)/:5+-(i-5p/^e*p)— ] < o 

lip pp p p 

(25) 
Tnk"^ 

{l+i*pPpy){TCT+Ti^T^d-{'^+i*TTpy){pCp+pi^pC^)+{l+i^pPpy) < 0. 

^ (26) 

Note that in absence of thermal conduction the above criteria become inde- 
pendent of the wave number k. 

Criterion (24) is Parker's ( 1953) isochoric criterion of thermal instability 
but modified by the presence of the chemical reaction, i.e. the change with 
time of temperature occurs in chemical equilibrium. 

Criterion (25) can be identified as Field's (1965) isentropic criterion for 
thermal instability also modified by the chemical reaction. In fact, as can be 
readily verified, if the equilibrium value .^*(p, T) is weekly dependent on den- 
sity and temperature, the relation (25) reduces to the corresponding Field's 
isentropic criterion. In the opposite case when the dependence of the heat- 
ing/cooling function on p and T is very week ( Cp ~ Ct ~ ), the criterion 
stands barely due to internal processes ( the chemical reaction ), i.e. 

T T Kk'^ 

(1 + TpiyU)[iM*p + — e*T)A + -(1 - Bpp^^p)—] < . (27) 

pp p p 

If {l+Tpu^^.T) > 0, the criterion simplifies further, that is, isentropic instabil- 
ity sets in when the expression in square brackets in relation (27) is less that 
zero. If the above occurs when (A^*p + ^C*r)>Cj < and (1 — Bpp^^p) > 0, 
sound waves with wave number smaller than the critical number kcdeHned 
by the equality in (27) are amplified and those larger ones are damped, but 
if {A^^p + -^^*t)j^^ < and (1 — Bpp^^p) < 0, sound waves are amplified 
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regardless of the wave length. On the contrary, if the expression in square 
brackets in (27) is positive, the sounds waves are damped. 

Relations (26) can be identified as Field's (1965) isobaric criteria for ther- 
mal instability, modified by the chemical reaction. As in the previous case, 
the respective thermal instability would proceed in chemical equilibrium. 
Generally, from the above modified criterion it follows that very fast chem- 
ical reactions may help to develop or may quench the respective isobaric 
thermal instability depending on the signs of ^^.t^^^ and ^^.pC^ terms for the 
particular reaction under consideration. In particular, if {l/fi + i^ppp) > 0, 
this criterion can be written as 

{TCt - PpCp) + {TU - (3pU)C^ + -y- < , (28) 

where 

^=#±1^. (29) 

Note that if the chemical parameter in chemical equilibrium weekly de- 
pends on density and temperature, i.e. ~ C,*p ~ 0, the above criterion 
reduces to the isobaric Field's criterion for nonreacting gases. However, the 
above generalized criterion, relation (28), provides further information. In 
effect, in the opposite case, that is, when the heating/cooling function jC 
weekly depends on p and T, Ct ~ >Cp ~ , criterion (28) simplifies to 

{Ti^T - (^pi.p)C^ + -y- < , (30) 

which can be identified as an isobaric instability criterion due only to the 
chemical reaction . 

Other important conclusion drawn from the instability criteria (24)- (26) is 
concerned with the critical wave numbers kc defined by the respective equality 
which corresponds to the marginal states. It is obvious that very fast chemical 
reactions may produce drastic changes in the critical scale lengths at which 
any thermal instability may occur. In particular, if Ct < 0, and ^^T-^g > 
but this last term is not large enough for quenching the isochoric instability, 
the critical wave number Ac ( = 27T/kc )can be considerable increased respect 
to the critical wave number in a non-reacting gas. Similar considerations can 
be made for the remaining two critical wave numbers, the isentropic and the 
isobaric ones, as it will be seen bellow where particular applications will be 
considered. 
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On the other hand, for fluctuations with wave number A; = 0, equation 
(20) becomes 

N^iaoN^ + aiN^ + + as) = , (31) 
where Oq and aistand as in equation (21), but 02 and 03 reduce to 

02 = (1 + -J^*tT) + Tc{2k[ + k^) ■ as = c{k[ + ks) . (32) 

If in addition t — , the dispersion relation (31 ) further reduces to 

N\a2N+a3) ^0 , (33) 

i.e. there is only one ( thermoreactive ) mode with rate A'' = —c{k[ + 
k'i)/{l + B^^^tT/A), which becomes unstable if k'-^ + /cs < 0, provided that 

1 + B^,tT/A > 0. 

For the opposite asymptotic case, i.e. when /c — > 00, the dispersion 
equation (20) simplifies to 

Af^ + -{2 + C*pPH-^ + ^(1 + C*pPH = . (34) 
T r 

Therefore, there are two modes with rates A/" = — 1/t and A/" = —(1 + 
i*pPP'^)/T, respectively. Obviously, for chemically stable reactions the last 
mode can be unstable provided that (1 + i^pPfii') < 0, otherwise any fluctu- 
ation with very short wave length is damped in a chemical time scale. 

3.2 General Case 

In the general case, when the time scale r is 7^ 0, case usually found in prac- 
tical problems, applying the Hurwitz criteria to the fifth order polynomial 
(20), five instability criteria are found. Explicitly, first instability criterion, 
Oi < 0, ( for T 7^ ) becomes 

(£,±£(M±s^J + !^!<0, (36) 

I I p 

with the sign -|- for chemical stable and the sign — for chemically unstable 
reacting fluids. The above criterion reduces to the isochoric criterion for 

nonreacting fluids when | r |— > cxd. On the other hand, for fluctuations 
with A; = 0, the instability criterion reduces to the condition [Ct ± R{2A + 
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^^,tBT)/ I r |] < 0, i.e. for a thermally unstable fluid ( £7- < ) the 
stabilizing effect of the chemical reaction increases when the chemical time 
scale I r | decreases, provided that 2 A + B^^tT > 0. 

Fifth instability criterion is just 05 < 0, and coincides with the isobaric 
modified criterion (28). Therefore, the above criterion is independent of the 
chemical time scale r and it stands as far as r 7^ 00, i.e. when chemical 
reactions proceed in the fluid under consideration. 

Second, third and fourth instability criteria become, respectively, 

A2 = aia2 - aoOa < , (36) 

A3 = A2a3 + (00^5 - 0104)01 < , (37) 

A4 = A3a4 + [— A2a2 + 00(0104 — ao05)]a5 < . (38) 
The relation (36) can be explicitly written as 

- ^)k' + ^[jT^*T + 2(2 + jT^.T)]}h 
+ -c^kl + -A;iA;3 + "-^k^ + ^{l + ^n,T)k^ < , (39) 

T T 7 T"' A 

which reduces to the isentropic Field's criterion when r — > 00. Therefore, 
it can be identified as a generalized isentropic criterion for thermochemical 
instability. For fluctuations with k — 0, but values of r 7^ 0, the above 
criterion simplifies to 



-(2+^re*T)(l+^re*T)+(4+3^re*T+2TcA;i)cA;i+(l+^Te.T+rcA;i)cA;3 < 0, 

(40) 

which in the asymptotic case when jCt ~ becomes 

(1 + jTi,T)[cks + i(2 + |re.T)] < , (41) 

or simply 

C.tA ± A(2^ + BTU) < , (42) 
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if (I + BT^^t/^) > 0, where + and — signs corresponds to chemically stable 
and instable reactions, respectively. 

Explicit relations for the criteria (37) and (38) are rather involved and 
generally, one must proceed numerically in specific applications. However, in 
the asymptotic case for k = 0, criterion (37) reduces to 

(2kl + 3 kiks + kl) c^kiT^+ 
\ A ^ A ^ J 

[^^'^ A ^ A A- + A- J^-° 

(43) 

If additionally ki ~ 0, the criterion simplifies to 

crk3[{l+ ^ )crk3 + {2+ — j — + — — — )] < (44) 

which can be identified as a thermoreactive criterion for instability. 

Fourth criterion (38) can also be identified as a second thermoreactive 
criterion for instability, as it can be verified considering the asymptotic case 
£p pa £t ~ 0, but £^ ^ 0. 

4 Astrophysical Applications 

4.1 CoUisionally Ionized Hydrogen Plasma 

This section will be devoted to applying the general results obtained in sec- 
tions 2 and 3 to the hydrogen plasma model studied in a previous work (IP), 
i.e. to a coUisionally ionized pure hydrogen plasma with a net rate function 
X{p, T,^) and a net cooling rate per unit mass jC{p,T,^), respectively, given 

by 

X{p, T, = NoaB{T)pe - Noq{T)pai - , (45) 
C{p, T, = NoRpeTPniT) + N^XP^^ - OiliT) + $(T)]-£o , (46) 
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where A^o is the Avogadro Number, x — 13.598 e.v, Cq is a constant per 
unit mass heating and the coefficients asiT)^ l^siT) and $(T) are 

given by Seaton (1959), Hummer and Seaton (1963) and Hummer (1963). 
The galactic value Cq = 3.25 x 10""^ ergs s^^ ( Potasch, Wesselius and 
Duinen 1979 ) has been taken as a reference value. Additionally, the thermal 
conduction coefficient is taken in the form 

«(p, T, = 2.5 X 10^(1 -e)TV2 + 1.84 xlO-^^-^^^^^ ^ (47) 

where the first term on the right hand side is the heat conduction by neutral 
atoms ( Parker 1953) and the second one is the Spitzer ( 1962 ) thermal 
conduction by electrons, where 



In A = 23.24 + liil(^)3(-^)i/2],j/r<4.2x 10=A', (48) 

'"^ = ''■'^ + '"' lQ»(ivIpO./2 l- '/I- > 4.2 X 10=if. 

The range of temperature under consideration is 3 x 10^ < T < 8 x lO^K . 

The equilibrium pressure as a function of particle number density has 
been plotted in Figure 1, on which the value of temperature of the marginal 
states have been also indicated. From solving the IP dispersion equation 
follows that regardless the value of k the hydrogen plasma is thermoreactiv- 
elly unstable in the range of temperature Ts ( = 1.66 x 10^ K )< T < T4 { 
— 4.30 X 10^ K ), with a maximum growing rate of 3.2 x 10~^ yr~^ at a tem- 
perature T = 3.2 X 10^ K . Additionally, there is an oscillatory mode which, 
for wave numbers A; < 2.1 x 10^^, is isobarically unstable in two temperature 
ranges , close to T ~ 9 x 10^ K and to T ~ 5 x 10*^ K , see Fig 2a. The lower 
range of temperature increases when k decreases. In Fig. 1 such a range Ti 
( = 8.65 xlO^ K )<T <T2 { = 1.26 x K ) is shown for the asymptotic 
value A; ^ 0. For A; > 1.7 x 10~^ this oscillatory mode stabilizes. The corre- 
sponding maximum growing rate is very short. For instance, a disturbance 
with k = lO^'', reaches a maximum growing rate -Re{nmax} = 7.8 x 10"^^ yr~^ 
a,tT — 9.66 X 10^ K at the lower temperature unstable range, and Re{n^sx} 
= 4.4 X 10~^^ yr~^ at T = 3.77 x 10^ K, for the parameter values used in 
Fig.l. Note that k is normalized to the value k^ = 1.89 x IO^Cq/RTiC^, where 
c* = (|i?Ti)^/^ and Ti = 157890i^r. Therefore, for the galactic value taken 
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by IP as a reference value, Cq — 3.25 x 10~^ erg one gets — IQ"^^ 

cm~^. For disturbances with larger wave numbers the oscillatory modes are 
damped. 

Figure 2b is a plot of the resulting rates from solving the dispersion 
equation (20). Qualitatively, the above results obtained by IP hold. However, 
the maximum values of the growing rates increase by a factor 2 and the fastest 
real damped mode ( for T > 10^ K , in Fig 2a ) becomes an oscillatory 
damped modc( for T > 3 x 10^ K , in Fig 2b ) . The above results also hold 
for larger wave numbers, as it can be seen in Figs. 2c and 2d , where the rates 
obtained from IP dispersion equation and from the fifth order polynomial 
(20), respectively, have been plotted as functions of k for a temperature 
T = 3.5 X 10^ K. The fourth order polynomial gives an oscillatory mode 
unstable for k < 2. x 10^'^{ critical value ) and two real modes, the slowest 
one being unstable. The fifth order polynomial also gives the above modes 
but increased by a factor 2, and an additional very fast damped real mode. 
The physical reason for the damping of the oscillatory mode at large enough 
values of k rests in the second viscosity originated by the chemical reaction 
( Landau & Lifshitz 1987 ). 

The effect of thermal conduction is twofold: it increases the damping and 
quenches the instabilities for any fluctuation with wave number k greater 
than a critical value kc- 

4.2 Photoionized Hydrogen Plasma with Metallicity Z 

In this section the results of section 3 will be applied to the photoionized 
hydrogen model studied by Corbelli and Ferrara ( CF ), i.e. an optically thin 

hydrogen plasma with metallicity Z heated and ionized by a background 
radiation field of mean photon energy E and ionization rate (. The net rate 
function X{p,T,^) and the net cooling rate per unit mass C{p,T,^) are 
respectively given by 

x(p,r,e) = Nopiea - (1 - e)e7c] - (1 - + 0)c, (49) 
£(p, T, = N^p[{i-c)zAHz+czAez+{i-m^eH+eKH+]-No{i-mEh+{i+m- 

(50) 

Except for the ionization rate indicated above by (, the remaining notation 
is like the one used by CF, as well as the corresponding expressions for: the 
number of secondary electrons (j) and heat released per photoionization 
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( ShuU & Van Steenberg 1985 ). the cooling efficiencies by collisions neutral 
hydrogen-ions and metal atoms K^z (Launay & Roueff 1977, Dalgarno and 
McCray 1972 ), electrons-ions and metal atoms Agz ( Dalgarno and McCray 
1972 ), Lya emission by neutral hydrogen Agij ( Spitzer 1978 ) and hydrogen 
recombination h-eH+i on the spot approximation (Seaton 1959 ). 

By solving numerically the fifth order dispersion equation (20), one ob- 
tains the following results: for a radiation field with a photon energy E = lb 
eV, a plasma with Z = 1 becomes stable in the whole range of temperature 
( T < 8.26 X lO^K ) where the thermochemical equilibrium may exist, con- 
trary to the CF result where there is instability in the range 33 < T < 703K. 
Figures 3a and 3b show the resulting rates, for a disturbance with k — 1 { 
in units of 10~^^cm~^ ), by solving the fourth [Eq. (2.7) of CF ] and fifth 
[Eq. (20) of Sec. 2 above ] order polynomials, respectively. In the first case, 
there is a damped mode ( the fastest one ) with two bifurcation points and 
an unstable oscillating mode in the above range of temperature. The disper- 
sion equation (20) leaves an additional non oscillating damped mode with 
the highest rate, unfolds the damped oscillating mode ( of Fig 3a ) in two 
damped real modes, and stabilizes the oscillating unstable mode (of Fig.3a 
) on the whole range of temperature. The rates corresponding to the above 
modes also are shown in Figs. 3c and 3d, as functions of the wave number 
k, for a temperature of T = 60K ( temperature at which becomes a 
maximum ). The slowest oscillating mode becomes unstable for k < 14.5 in 
Fig 3c; instead, it remains stable in Fig. 3d with a rate being an increasing 
function of k. The maximum damping occurs for acoustic oscillations with 
k > 10^. This effect is expected to occur as a consequence of the second 
viscosity due to the chemical reaction (Landau and Lifshitz 1987). Addition- 
ally, the new damped non- oscillating mode shows dispersion at high values 
of /c ( > 3 X 10^). Similar behavior is found at other values of temperature. 

Fig 4a (fourth order polynomial ) and 4b (fifth order polynomial ) show 
the rates as function of temperature for disturbances with A; = 1, in a plasma 
with Z — 1, and photon energy E — 10^ eV. With the exception of appear- 
ing the fifth non oscillating damped mode ( with a rate value between the 
corresponding rates values of the two modes given by the fourth order poly- 
nomial ), no qualitative changes appear. Only slight changes are observed in 
the values of temperature at which bifurcation occurs for both modes, the 
fastest stable and the slowest unstable one. In Figs 4c (fourth order polyno- 
mial ) and 4d (fifth order polynomial ) the corresponding rates are plotted as 
functions of the wave number k, for a temperature T = 200i^. The rates of 
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the damped mode and the corresponding value of k where bifurcation occurs 
are shifted to high vahics while the unstable oscillating mode ( Fig 3c ) is 
stabilized for wave numbers k > 4.8 ( Fig.4d ). 

Figs. 5a, 5b are as Figs. 3a, 3b but for Z = 0.5; and Figs. 5c, 5d are as 
Figs. 3c and 3d but for a temperature T — 500K . There, in the whole range 
of k under consideration the fourth order polynomial gives two oscillating 
modes, the fastest one being stable and the slowest one unstable for k < 4.8 
( Fig. 5c ). Instead, the fifth order polynomial unfolds the above damped 
oscillating mode in two non oscillating modes and stabilizes the unstable ( 
slowest ) mode. Additionally, the new non- oscillating damped mode also 
shows dispersion at high wave values of A; ( > 10^ ), and the acoustic mode 
is also strongly damped by second viscosity for values of /c > 2 x 10^ ( Fig.5d 
). Similar results hold for other values of the temperature. 

For low metallicity plasmas ( say Z = 10^^ ), the fifth order polynomial 
gives three non oscillating and one oscillating mode, or two oscillating and 
one real mode. In particular, for E — lOOeV, and in the range temperature 
where the plasma is unstable, i.e. 72 < T < 4.34 x lO^K, the fourth order 
polynomial gives a stable and an unstable real mode ( the slowest one ), and a 
damped acoustic mode. Instead, the fifth order polynomial gives an unstable 
real mode ( also the slowest one ) and two damped acoustic modes. This 
unstable thermal mode is slightly faster than the corresponding one resulting 
from solving the fourth order polynomial. On the other hand, for Z — 0, the 
plasma also becomes thermochemically stable. 

Specifically, in clouds located outside the optical disk of galaxies or in 
the galactic halo the probable values for the energy and metallicity, quoted 
out by CF are: = 20 , Z = 1 and T 800 K. For these values the 
corresponding rates have been plotted in Figures 6a-6b ( as functions of T ) 
and 6c-6d ( as functions of A; ). The range of temperature for which oscillatory 
instability may set in according to the CF analysis is rather wide 43 < T < 
3.21 X 10^ K, ( Fig. 6a ) , instead, the present analysis gives a considerable 
reduction of it, i.e. the plasma is unstable only in the range of temperature 
49 < T < 8.46 X 10^ K. Additionally, while the fourth order dispersion 
equation gives oscillatory instability for any value of the wave number and 
maximum growing rates ~ 10~^ yr^^ , for disturbances with k > 10~^^cm~^( 
Fig 6c ), from the dispersion relation (20) one obtains that there is instabihty 
only for fluctuation with k < kc { Fig. 6d ), where the critical value kc is 
weakly dependent on temperature ( 10~^^ < kc < 9.6 x 10~^* cm~^ ), in the 
above range of temperature where the plasma is unstable ). The maximum 
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growing rate is strongly dependent on T, and the upper value rimax = 4.8 x 
10~^?/r~^is reached at temperature T ^ 500 K for fluctuations with wave 
number k = 6.0 x 10^^^ cm~^ (see Fig.6d , dashed line ). 

^From the above , one may conclude that the effect of taking into ac- 
count the change of the ionization equilibrium position due to the variation 
of density and temperature, is to stabilize the oscillating modes, at least for 
high enough values of the wave number, as it should be expected from Lan- 
dau's results on second viscosity. However, the quantitative stabilizing effect 
depends on the particular values of photon energy , metallicity and plasma 
temperature. For instance, for the values used in Figures 3c-3d and 5c-5d, 
such an effect is able to stabihze the unstable acoustic mode in the whole 
range of k under consideration. But, for those used in Fig4c-4d the effect is 
only capable of stabilizing acoustic oscillations with k > 4.8. 

5 Summary and Conclusions 

In summary, if one takes into account the effects of a finite chemical relaxation 
time in optically thin reacting plasmas a fifth order dispersion equation is 
obtained, instead of the fourth order polynomial obtained when such effects 
are neglected. The quantitative effects on the previous four modes depend 
on the particular reacting plasma under consideration, but generally, it tends 
to further stabilize the plasma and to increase the damping of stable modes 
because of the second viscosity generated by the chemical reaction. 

In the coUisionally ionized pure hydrogen plasma model previously stud- 
ied by IP the resulting effects are not too strong. In addition to the ap- 
pearance of the new damped real mode, only an increase in the value of 
both the damping and growing rates occur. However, in the photoionized 
plasma studied by CF, the effects may result more severe, depending on the 
exact value of the parameters under consideration. In particular, the range of 
temperature where the plasma is unstable can be strongly reduced, and the 
oscillatory instability can be quenched for k large enough. So for instance, 
for values likely occurring in the galactic halo, the present analysis leaves a 
very restrict values of k for developing instability. See the end of the previous 
Section. 
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FIGURE CAPTIONS 

Fig. 1. The equilibrium pressure as a function of particle number den- 
sity for a collisionally ionized pure hydrogen plasma. The temperatures of 
marginal states are indicated by the labels T1...T4, see text. 

Fig. 2. The resulting rates from solving the fourth order polynomial of IP 
as functions of temperature for a perturbation with wave number k = 10~^ 

(a) , and as functions of k for a temperature T = 3.5 x 10*^ K (c). The resulting 
rates from the fifth order polynomial (20) as functions of T for k = 10^^ (b), 
and as functions of A; for T = 3.5 x 10"^ K (d). For a heating rate value of 
the order of the galactic heating rate, k is in units of 10~^^ cm~^. Positive 
rates of complex modes are indicated by C"*" and those corresponding to real 
roots by R^. 

Fig. 3. The rates resulting from solving the fourth order (a) an fifth order 

(b) polynomials ( see text ) as functions of temperature for, A^oP = 1 cm^^, 
a photon energy £■ = 15 eV, metaUicity Z — 1, and k — 1. The wave number 
k is in units of 10~^^ cm~^. The values of the two negative intermediate 
real rates are indistinguishable at the scale of Fig. (b) . Figures (c) and (d) 
respectively correspond to the rates on Figures (a) and (b) but as functions 
of k, and for T = 60 K. As in Figure (b) the two damped intermediate real 
roots are also indistinguishable at the scale of Fig. (d). The range where the 
oscillating mode is unstable is indicated by C+. 

Fig. 4. Figures (a) and (b) are as Figures 3a and 3b, respectively, but for 
E = 100 eV . Figures (c) and (d) are as Figures 3c and 3d, respectively, but 
for E = 100 eV and T = 200 K. The range where the real mode is unstable 
is indicated by R^. 

Fig. 5. Figures (a) and (b) are as Figures 3a and 3b, respectively, but 
for Z = 0.5 . Figures (c) and (d) are as Figures 3c and 3d, respectively, but 
for Z = 0.5 and T = 500 K. 

Fig. 6. Figures (a) and (b) arc as Figures 3a and 3b, respectively, but 
for £■ = 20 eV, Z — 1 . Figures (c) and (d) are as Figures 3c and 3d, 
respectively, but for E ^ 20 eV, Z ^ 1 and T ^ 800 K. On Fig (d) the 
dashed curve corresponds to a temperature T — 500 K, value at which the 
maximum growing rate is reached. 
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